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CONSIDEWd^IOWG FOK T1{E APPI.ICATION OF FINITE-ELEf.ffiKT BEAM MODELING 
TO VIBRATION ANALYSIS OP FLIGHT VEHICLE STRUCTURES* 

Raymond G. Kvaternik 
Langley Research Center 

SUNKARY 

The manner of representing a flight vehible structure as an assemhly of 
heam, spring, and rigid-body components for vibration a:ialysis is described. 

The development is couched in terms of a substructures methodology ■which is 
based on the finite-element stiffness method. The particular manner of employ- 
ing beam, spring, and rigid-body components to model such items as wing struc- 
tures, external stores, pylons supporting engines or external stores, and 
sprung masses associated with launch vehicle fuel slosh is described by means 
of several simple qualitative examples. A detailed numerical example consist- 
ing of a tilt-rotor VTOL aircraft is included to provide a unified illustration 
of the procedure for representing a structure as an equivalent system of beams, 
springs, and rigid bodies, the manner of forming the substructure mass and 
stiffness matrices, and the mechanics of writing the equations of constraint 
which enforce deflection compatibility at the j-unctions of the substructures. 
Since many structures, or selected components of structures, can be represented 
in this manner for vibration analysis, the modeling concepts described and their 
application in the numerical example shown should prove generally useful to the 
dynamicist. • 

*The information presented herein is based on a portion of a thesis 
entitled "Studies in Tilt-Rotor VTOL Aircraft Aeroelasticity," which was sub- 
mitted to Case Western Reserve University, Cleveland, Ohio, in partial fulfill- 
ment of the requirements for the degree pf Doctor in Philosophy in Engineering 
Mechanics, J'une 1973. 
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INTRODUCTION 


Customary engineering practice in dynamic analyses of complex aerospace 
structures consists of introducing some type of physical idealization to 
simplify the structure, establishing a finite-element mathematical model of 
the idealized system, and solving the resultant equations of motion. Finite- 
element models required for static analyses are usually quite complex in order 
to predict the stresses satisfactorily. For dynamic analyses, however, simpler 
models based on the use of (one-dimensional) beam elements to represent either 
the entire structure or selected components are often adequate (refs. 1 to ll). 

A self-contained treatment of the practical aspects of applying finite- 
element beam modeling techniques in dynamic analyses is, for the most part, 
either lacking in the literature or of a fragmentary nature. The purpose of 
this report is to present a unified treatment of several aspects of finite- 
element beam modeling as might be applied to vibration analyses of flight 
vehicle structures. Herein, beam modeling is used in the general sense to 
include the use of beam, spring, and rigid-body components. The discussion is 
couched in terms of the substructures approach to vibration analysis as devel- 
oped in reference 12. Several simple qualitative exajnples illustrating the 
manner of employing beams, springs, and rigid bodies to model such items as 
wing structures, external stores, pylons supporting engines or stores, and 
sprung masses associated with launch vehicle fuel slosh are given. A detailed 
numerical example consisting of a tilt-rotor VTOL aircraft structure is employed 
to illustrate the manner of representing a structure as an equivalent system of 
beam, spring, and rigid-body substructures, the formation of the substructure 
mass and stiffness matrices, and the mechanics of writing the equations of 


2 


constraxnt which enforce deflection corapatihility at the junctions of the sub- 
structures. The direct method of vibration analysis as described in refer- 
ence 12 and implemented in a special-purpose computer program designated SUDAN 
(substructuring in Direct ANalyses) is used to solve for the modes and 
frequencies. 

The present report relies on reference 12 for a complete understanding of 
the theoretical basis of the procedures implemented in the SUDAN program. For 

this reason it is recommended that reference 12 be used in conjunction with 
this report. ’ 




SYMBOLS 


[A],[B] 

AE. 

1 

[C] 


e. 

1 


El. 

1 

GJ. 

1 

I. 

1 

I. . 

IJ 

[K] 

[Klfi 

^^^B-S 


[AK] 

[K]^ 

im 


[k] 


(i) 


K .K, 


M 

[M] 


submatrices of 

extensional rigidity of ith beam element 

submatrix of [K]gj also matrix of coefficients of constraint 
equations expressed in discrete physical coordinates 

offset of ith mass from beam elastic axis 

flexural rigidity of ith beam element 

torsional rigidity of ith beam element 

lumped torsional inertia at ith beam station 

components of rigid body inertia tensor (See eq. 15* ) 

composite matrix containing free-body stiffness matrices of 
substructures as submatrices on the principal diagonal 

stiffness matrix for beam bending 
stiffness matrix for beam-spring 

stiffness matrix for beam torsion 

matrix of additional spring stiffness terms 
stiffness matrix for beam extension 

com.posite matrix containing the free— body modal stiffness 
matrices of substructures as submatrices on the principal 
diagonal 

stiffness matrix of ith substructure regarded as a free-body 
longitudinal and lateral hub flapping spring rates 

lengths; length of ith beam element 
mass of rigid body component 

composite matrix containing free-body mass matrices of substruc 
tures as submatrices on the principal diagonal 

1+ 




lumped mass at ith beam station 

[mJb 

mass matrix for beam bending 

[Mlj 

mass matrix for beam torsioij 

[M] 

III 

mass matrix for beam extension 

[Jl] 

composite matrix containing the free-body modal mass matrices of 
substructures as subraatrices on the principal diagonal 

^'^^b+rb 

beam + rigid body mass matrix 

[mIb+t 

coupled bending— torsion mass matrix 

[M]c_,p 

mass matrix of wing carry-thrqugh structure 

ROTOR 

mass matrix of rotor 

PYLON 

mass matrix of pylon 


mass matrix of rigid body 


mass matrix of ith substructure regarded as a free body 


column matrix of discrete, physical coordinates for ith sub- 
structure regarded as a free body 

q 

generalized coordinate 

{q},{^^ 

jth nodal degree of freedom 

colimm matriic containing the discrete, physical coordinates of 
all the substructures 

[R],[T],[L] 

connection coordinate transformation matrices 

RI. 

X 

RI 

lumped rotary inertia at ith beam station 

rotary inertia of rigid body about its center of gravity 

T 

kinetic energy 
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kinetic energy of ith wing section 


x,y,z 


X ,y ,z 
r»*'r r 


a, 3 ,Y 


a ,3 ,Y 
r r 'r 


iincoupled system modal expansion matrix containing substructure 
modal subsets on principal diagonal 

extensional displacement of beam 

extensional displacement of beam at ith station 

potential energy 

bending displacement of beam 

bending displacement of beam at ith station 

column vector of jimction coordinates, r = 1, 2, . . . , 15 

column vector of auxiliary coordinates, s = i, ii, , , vii 
displacements 

connection coordinate displacements 
rotations 

connection coordinate rotations 
proprotor disc coordinate vector 
bending slope of beam 

pylon tilt angle forward from vertical position 

angles; bending slope of beam at ith station 

column matrix containing all substructure modal coordinates 
twist of beam at ith station 


Physical quantities in this report are given in the International System of 
Units (si). The calculations were made in Customary Units. 




SUBSTRUCTURES APPROACH TO DYNAMIC ANALYSIS 

Three distinct steps may be identified with dynamic analyses based on a 
substructiires approach: partitioning into substructures, discretization of the 

substructures, and assembly of the mathematical models of the substructures. 
Structural partitioning consists of dividing the structure into a collection of 
smaller components or substructures which can be more easily handled. An illus- 
tration of such a partitioning is depicted in Figure 1. Discretization involves 
establishing a finite-element mathematical model of each isolated substructure 
based on the known inertial and elastic properties of the substructures. Detailed 
considerations relating to finite-element modeling are available in several 
books. (See, for example, references 13 to 17.) Application of the finite- 
element stiffness matrix method to each substructure regarded as a free body 
leads to a discrete mass matrix [m], a discrete stiffness matrix [k], and a' 
vector of discrete coordinates {p} for each substructure. The free vibration 
equations of motion corresponding to each substructure can be collected into one 
partitioned diagonal matrix equation having the uncoupled form 




, 

r •> 

Cp}'^ 




< 



) 

___ 



^ J 





J 


(0)(NS) 


^Discrete coordinates define the translations and rotations at a set of 
discrete points. 
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where US is equal to the number of substructures into which the structure has 
been divided. For notational convenience equation 1 is written in the compact 
form 

[M] {z} + [ic] {z} = {0} (2) 

For the substructures shown in Figure 1, i’or example, the composite matrices 
[M] and [K] appearing in equation 2 each would have the form given diagrara- 
matically in Figure 2. Each block in Figure 2 represents the mass or stiffness 
matrix appropriate to a substructure. The ordering of the submatrices within 
[M] and [K], though required to be consistent, is arbitrary. For component 
mode synthesis the modal equivalent of equation 2 is 

[JI] ih + M] iO = {0} (3) 

where [^] and [xK] are the composite matrices containing the substructure 
modal mass and stiffness matrices as submatrices on the principal diagonal. 

A consequence of any substructuring procedure is the introduction of coor- 
dinates which are not generalized coordinates but are related by equations of 
constraint which must be imposed to restore geometric compatibility at the 
interfaces of the substructures. Since the matrices [M] and [K] in equa- 
tion 2 and and [‘3<] in equation 3 are established on the basis of such 

a substructuring procedure, neither the coordinates forming the vector {z} nor 
those forming the vector {C) are independent. The dependency equations relat- 
ing the various discrete coordinates in {z} can be put into the matrix form 

[C] {z} = {0} ik) 
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and those relating the various modal coordinates in {?) ean be put into the 

form 


[C] [U] iO = {0} 


( 5 ) 


Vhare fc] is a constant matrix depending solely on the geometric configuration 

Of the interfaces and [U] - the uncoupled system modal expansion matrix. A 

set Of independent coordinates consistent vlth the appropriate equations of con- 
straint must be established. Usual practice vhen dealing with equations of 
constraint (see ref. 16. for example) is to select certain of the coordinates 
as independent coordinates and express the remaining dependent coordinates in 
terms of those which have been selected to be independent by solving the con- 
straint equations as simultaneous equations. A different approach to establish- 
ing independent coordinates in the presence of equations of constraint is an 
algorithm (ref. 19) which is based on solving an eigenvalue problem associated 
with a symmetric matrix formed from the coefficients of the constraint equations. 
The basis of these procedures is summarized in reference 12. 

BEAM MODELING OP SUBSTRUCTURES 
General Considerations 

Since the substructures are treated as distinct and separate components in 
a substructuring methodology, their structural properties are most conveniently 
defined relative to axes local to each component. The specification of the 
mass and stiffness matrices corresponding to each of the three types of struc- 
tural members employed in the be™ modeling of a structure is the subject of 
the following subsections. 
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Su~b structures ~breated as ~beaniSt— The elastic properties of substructures 


treated as beams are teiken to be defined in terras of the distribution of 
flexural, torsional, and extensional stiffness (EJ, GJ , and AE, respectively) 
along a theoretically defined elastic axis. The continuous distortion of the 
beam in extension, torsion, and bending in two perpendicular planes is approxi- 
mated by specifying appropriate deflections and rotations at a number of dis- 
crete points or stations along the elastic axis of the beam. A beam segment or 
element is defined to be the length between two such stations, The stiffness 
of each element is assumed to be constant and given by the average value of the 
stiffnesses at the two adjacent stations. The distributed mass of the beam is 
discretized by replacing the distributed mass within mean locations on either 
side of the stations established by statically equivalent concentrated masses 
at each of the stations along the beam. In general, each mass has three com- 
ponents of translational inertia and three components of rotational inertia. 

Mass static unbalance about the local elastic axis is also preserved. 

The mass and stiffness matrices of an unrestrained, arbitrarily oriented 
beam element are each of order 12 12, three translational and three rotational 

degrees of freedom being associated with each of its two ends. If the local 
coordinate axes are chosen to coincide with the principal axes of the cross 
section, the 12 x 12 stiffness matrix can be expressed in terms of uncoupled 
4 X li and 2 x 2 submatrices located on the principal diagonal (ref. 13). A 
similar partitioning is possible for the mass matrix if the centers of gravity 
of the sectional masses are located on the elastic axis of the beam. In general 
however, the sectional centers of gravity of the original component will not lie 
along its elastic axis. This will lead to additional, off-diagonal, mass terms 
which will couple bending and torsion and/or bending and extensional motion. 
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It will 1)6 assumed tha,t the sectional centers of gravity lie on the elastic 
axis for the present and later in the report it will be shown how to accotmt 
for any mass coupling terms, The form of these submatrices for a beam sub- 
structure is shown below fop the two-element beam shown in figure 3. 

Beam bending .- The stiffness matrix for vertical bending is put into the 
partitioned form: 



where, for a beam modeled using two elements, 



11 




The hending displaGeiaentis and slopes have each been grouped together and placed 
in the order shovm in equation 6 for convenience in their computer implementation. 
Extension to additional elements is apparent in the distinctive forms of [a], 

[B] , and [C]. The corresponding inertia matrix is given by 





where the matrix elements not shown are zero. Similar matrix expressions des- 
crihe lateral bending. 


Beam torsion: The torsional stiffness matrix is given by 



and the torsional inertia matrix by 


?-3 





( 12 ) 



H 

1 


♦i 

Mt = 

^2 






*3 


Extension to piore elements is olDvious. 

Beam extension: The stiffness matrix describing extensional (axial) 

deformation has the form 


M.- 


AE, 


AE, 


L 


AE, 




Ail. 


AE, 


AEg 


J 


(13) 


and the corresponding inertia matrix is 

r 


Q- 


^2 ^3 


[m] = 


A 


^1 




2 


2 



^3 


(lit) 


Again, the extension to additional elements is obvious. 





Su'bstructures treated as springs .- In many instances actual springs may 
comprise some of the components of a structure or members may be treated as 
one- or two«-degree-of- freedom springs depending on whether they have one end 
tied to ground or both ends free. These springs are defined by 1 x 1 or 2 x 2 
stiffness matrices having spring constants as matrix elements and associated 
1 X 1 or 2 X 2 inertia matrices which a,re null (zero). 

When modeling short structural components as beams it is sometimes con- 
venient to account for their elastic characteristics in the definition of the 
substructure but to include their inertial properties with adjacent components, 

A massless, uniform beam segment is employed for this purpose. Since only the 
elastic properties of this beam segment are considered it has the character of 
a spring and it seems appropriate to include its description here. For descrip- 
tive purposes this massless beam segment will be referred to as a "beam-spring" 
herein. The 12 x 12 stiffness matrix for the beam-spring having the coordinate 
ordering shown in figure U has the form given in figure 5. The corresponding 
inertia matrix is taken to be null. 

Substructures treated as rigid bodies .- Components such as ordnance, 
external fuel tanks, engine/nacelle combinations, etc., can often bo treated as 
rigid in dynamic analyses. The inertia matrices of such rigid bodies with 
respect to local body axes at the center of gravity have the general form 
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where the elements not shown are zero. The corresponding stiffness matrices 
are null since a rigid body has no strain energy associated with its motion. 

Oftentimes, while it may be necessary to account for the translational or 
rotational rigid-body motion of a beam substructure in some direction, the cor- 
responding deformation can be neglected. This situation is easily accommodated 
by this rigid-body component. 

Once the mass and stiffness matrices for the substructures have boen 
determined, the mass and stiffness matrices for the complete structure are given 
by the composite matrices formed by locating the substructure matrices along the 
principal diagonal. These composite matrices are denoted [M] and [k] , 
respectively. For the "stick" model shown in figure 6, for example, [m] and 
[K] would have the form shown in figure 7. It should be noted that the stick 
model of figure 6 and the particular freedoms indicated in figure 7 are appropri- 
ate to a symmetric vibration analysis. Each block in figure 7 corresponds to a 
submatrix. The ordering of these submatrices within the larger substructure 
submatrices (indicated by braces) and the ordering of substructure submatrices 
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within [M] and [K] must be compatible but is otherwise arbitrary. Since the 
mass and stiffness matrices for each substructure are generated independently, 
no inter- substruct-ure coupling exists in [m] or [K]. However, intra- 
substructure coupling (i.e., coupling between subn^atrices within a substructure 
submatrix through off-diagonal terms) can exist. Situations in which such 
coupling arises will now be discussed. 

Specific Considerations 

If the sectional centers of gravity of the lumped masses for the beam sub- 
structures do not fall on the elastic axis, there will arise mass static unbal- 
ance terms which will couple two or more of the blocks in the mass matrix for 
the beam substructure. For example, if the sectional centers of gravity of the 
wing or tail surfaces of the airplane in figure 6 were displaced from their 
elastic axes in the plane of the surface the mass static unbalance terms would 
appear outside the block diagonal areas indicated in figure T and couple the 
vertical bending and torsion submatrices in the matrices fop the wing or tail 
substructures in [M]. Mass coupling terms of this type will also occur if a 
rigid body component is not treated as a unique substructure (as discussed 
earlier) but has its mass matrix combined with the mass matrix of the beam com- 
ponent to which it is attached. If, for example, the main landing gear assembly 
of the aircraft shown in figure 6 were treated as a rigid mass and its inertial 
properties combined with the forward fuselage beam, the vertical bending and 
axial rigid body submatrices of the forward fuselage beam substructure 
in [M] would be coupled.* Coupling terms would also arise within the vertical 


*For an anti-symmetric analysis the lateral bending and torsion blocks of 
the substructure would be coupled. 
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bending block, leading to a non-diagonal mass matrix for beam bending. Inter- 
and intra- substructure coupling terms arise in the stjLffness matrix of the 
partitioned structure, [K], when the stiffness characteristics of springs 
which are not treated as substructures are combined with the stiffness matrices 
of members to which they are attached rather than treated as separate substruc- 
tures. These aspects, as well as pthers , are elucidated below in several 
simple qualitative examples. 

Treatment of ma ss static unbalance .- Consider the sectionalized wing plan- 
form shown in figure 8. The section (lumped) masses, bending inertias, and 
torsional inertias (about the c.g.) are denoted by M^, and I., respec- 

tively. The perpendicular distance between the seption c.g. locations (assumed 
to be in the wing chord plane) and the wing elastic axis are denoted by e.. 

The kinetic energy of each section, T^ , expressed in terms of displacements 
and rotations of the elastic axis station, has the form indicated at the bottom 
of figure 8, Substituting this expression into Lagi'ange’s equation 

d^/BTN 9V _ ^ 

" 9q ~ ° ( 16 ) 

and performing the appropriate differentiations leans to the coupled, bending- 
torsion mass matrix shown in figure 9. Assuming figure 8 is appropriate to the 
wing of the aircraft in figure 6, then [m] has the form shown in figure J. 

The first two blocks of substructure would be given by twice the mass matrix 

of figure 9- If the sectional center of gravities were also displaced vertically 
from the wing elastic axis (i.e., not in the wing chord plane) additional coupl- 
ing terms would arise and the vertical bending, torsion, and fore-and-aft 
bending blocks of substructure ffk would be coupled. 
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Alternative treatment of rigid bodies .- As an alternative to treating a 
rigid body as a substructure, the inertial properties of the rigid body (total 
mass lumped at the c.g. and moments and products of inertia relative to axes 
fixed in the body at the c.g.) can be combined with the inertia matrix of the 
elastic component to vhich it is attached. This treatment has the advantage of 
not explicitly introducing the degrees of freedom associated vith the rigid 
body into the problem. 

Rigid attachment: If the rigid body is rigidly attached to an elastic sub-r 

structure, the procedure consists in modifying the kinetic energy expression for 
the substructure to include the effects of the concentrated mass, inertia, and 
static unbalance about its point of attachment. For example, consider the 
situation depicted in figure 10. Suppose one wishes to combine the inertial 
properties of the center-^line fuel tank, regarded as a rigid body, with the 
inertia matrix of the fuselage beam by regarding the tank as rigidly connected 
to the nth fuselage beam station. Assume that each fuselage mass has both 
vertical bending and longitudinal degrees of freedom. If the principal inertia 
axes of the fuel tank are parallel to the principal geometric axes of the fuse- 
lage beam the kinetic energy of the tank expressed in terms of the motion at 
the nth mass station on the beam has the form given at the bottom of figure 10. 
Substituting this expression into Lagrange's equation gives the matrix of addi- 
tional tezms which must be added to the mass matrix (assumed to be diagonal) for 
the fuselage beam. The final mass matrix is given in figure 11. 

Flexible attachment: If the fuel tank of figure 10 were attached to the 

fuselage beam through a flexible member which could be treated as a spring sub- 
structure the inertia properties of the tank could be combined with the (null) 
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inertia matrix of the spring substructure using a procedure similar to that 


described directly above. 

Alternative treatment of springs .-. For convenience, springs can be divided 
into those which have one end tied to ground (one-degree-of-freedom springs) 
and those which have both ends ’’free" (two-degree-of-freedom springs) in the 
sense that while both ends are attached to some component, neither end is tied 

to ground. 

Springs having one end tied to ground: If the spring attachment point on 

the structure is on the clastic axis at a point identified as a degree of free- 
dom, the spring constant(s) can be simply added to the appropriate diagonal 
term(s) of the stiffness matrix of the component. If the point of attachment is 
on the elastic axis but not at a point Identified as a degree of freedom, coupling 
»111 occur, in the latter case It Is sometimes convenient to introduce an aux- 
iliary massless station at the spring attachment point. Similar considerations 
apply to the case where the spring is connected to the elastic axis through a 
rigid offset. These situations ate depicted in figure l?a. The general form of 
the additional stiffness terms due to springs attached to the elastic axis at a 


station identified as a degree of freedom is indicated in figure i3a. 

Free-free springs; The spring constants of springs which have neither ana 

tied to ground, such as given in figure 12h, can he combined with the stiffness 
matrices of the substructures to which they are attached by writing the potential 
energy of the springs in terms of coordinates at the points of attachment. The 
general form of the additional stiffness terms due to the spring coupling of 
figure 12b is shown in figure 13b. A specific illustration of the use of this 
expedient in the realm of launch vehicle dynamics may he given with the aid of 

figure 14. 
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In dynamic analyses of launch vehicles the dynamic effects of sloshing 
propellants are usually included by introducing a dynamically equivalent mechani 
cal analogy , composed of fixed and oscillating masses connected to the tank by 
springs or pendulums, to account for each important vibration mode of the liquid 
as a degree of freedom (ref. 20). This equivalent lumped-parameter mathematical 
model can then be combined with appropriate discrete-element representations for 
onher components of the vehicle. A spring-mass analogy is shown in figure l4. 
One such spring-mass assembly is provided to represent the dynamic effects of 
sloshing accompanying vertical bending (translation and rotation) and longitudi- 
nal oscillation. For illustrative simplicity, all three sloshing masses are 
taken to be attached to the same (nth) beam station. The non-sloshing portion 
of the fluid would simply be combined with the beam mass at the nth station 
of the beam. The potential energy of the springs, expressed in terras of the 
deflections of the nth beam station (w^, 0^, and the deflections of the 

slosh masses (w, 0, u), is given by V in the figure. Substituting this 

expression into Lagrange's equation leads to the matrix [AK] shown in fig- 
ure 15. [AK] is the matrix of spring stiffness terms which must be added to 
the stiffness matrix for the beam. 

ILLUSTRATIVE EXAMPLE 

An example illustrating the manner of modeling an aircraft structure using 
beam, spring, and rigid body components, the manner of forming the substructure 
mass and stiffness matrices, and the mechanics of writing the equations of con- 
straint for a structural configuration of some engineering complexity is 
presented in this section. The subject configuration is the tilt-rotor VTOL 
aircraft design shown in figure l6,: It should be emphasized that the 
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ideallizatiQn of a complete aircraft as a simple "stick" model using a beam-type 
representation for all the structural components should be applied with engi- 
neering judgement. Such a simplification may be inadequate to account for the 
influence of such things as large structural cutouts in the airframe (e.g., 
doorways) and highly redundant structural components (e.g., wing carry-tiirough 
box). The intent of the following example is solely to provide an illustration 
of the practical aspects of the application of beam modeling which can be 
readily adapted to other configurations. 

Beam Model Eepresentation of Aircraft 

The subject aircraft is depicted in silhouette form in figure IT along with 
the idealized model established using the beam, spring, and rigid-body compo- 
nents described above. The fuselage, wing, and empennage structures are 
replaced by non-uniform beams lying along the theoretical elastic axes of the 
respective components. Since the fuselage elastic axis has two changes in slope, 
three beams are used to represent the fuselage structure. The root ends of the 
wing and vertical tail beams are located at the pei*iphery of fuselage cross- 
section. The wing carry-through structure joining the root ends of the wings 
is idealized as a heam-spring, its (rigid-body) inertial properties being com- 
bined with the inertia matrix of the second fuselage beam. The geometric offset 
between the second fuselage beam and the wing carry— through beam— spring repre- 
sents the fuselage depth between the, fuselage elastic axis and the carry-through 
elastic axis. The transition fi’om the fuselage beam to the wing carry-through 
beam is made by equations of constraint to express a rigid connection. A 
similar treatment is employed for the geometric offset between the aft fuselage 
beam and the vertical tail beam. The pylon otru.cture, consisting of the 
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transmission/ engine assembly, is treated as a rigid body while the mast (the 
portion of the drive shaft extending forward from the transmission) on which 
each proprotor is mounted is treated as a beam-spring. The geometric offset 
between the wing elastic axis at the wing tip and the position on the trans- 
mission centerline through which the conversion axis about which the pylon tilts 
passes is specified to be rigid via equations of constraint. The proprotor 
blades are assumed to be rigid. Since the blades are rigidly attached to the 
hub in a gimbaled proprotor design (such as the design depicted in fig. l6) the 
proprotors (assumed to be non-rotating) are treated as rigid discs in the 
analysis. The hub is initially taken to be rigidly attached to the mast. How- 
ever, for illustrative purposes, the manner of treating a hub which is spring- 
connected to the mast will also be described. 

Since the substructures are treated as independent components, their 
structural properties and deflection characteristics are most conveniently 
defined relative to axes local to each component. The existence of such local 
axes for the definition of these quantities is assumed herein. The local com- 
ponent axis systems are also employed to establish the deflection compatibility 
equations at the junctions of the substructures. The set of right-handed local 
junction coordinate axes used in writing the equations of constraint is identi- 
fiable by subscripts in figure 18. Each of the directions so indicated is taken 
to be positive. Vectors representing positive rotations (not 

shown) about x^, y^, z^, respectively, are taken in the same direction as 

vectors representing positive x^, y^, z^. A shorthand notation for a column 

vector of these junction, or connection, coordinates is given by where 
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Several auxiliary right-handed coordinate systems (x'fhich are not associated with 
degrees ol freedom) are employed to facilitate writing the constraint eoLuations. 
These auxiliary coordinate systems are also shown in figure l8 and are distin- 
guished by Boman numeral superscripts. A shorthand notation for a column vector 
of these auxiliary coordinates is given by {x} where 



The structural data for the fuselage, wing, and empennage of this design 
were available (x“ef. 21.) in the form of bar graphs showing the distribution of 
mass, rotary inertia, and mass static unbalance and curves showing the variation 
of El and GJ along the calculated elastic axis of each member. The wing and 
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tail inertial data were "based on cuts perpendicular to the elastic axis. Data 
for the proprotor and pylon consisted of lumped mass properties and selected 
EI/GJ characteristics. Employing these data, lumped-mass/ stiff ness discretiza- 
tions were established for the fuselage, wing, and empennage beams. These are 
summarized in Tables 3. to 3. Inertial properties of components treated as rigid 
bodies and beam-springs are given in Table 4 while the elastic properties of 
beam-spring components are listed in Table 5. A summary of the pertinent geo- 
metric quantities is contained in Table 6. Utilizing aircraft symmetry about a 
vertical plane through the center of the fuselage, explicit consideration need 
be given to only one-half of the aircraft. A consequence of this separate 
treatment of the symmetric and anti-symmetric problems is the identification of 
displacements and rotations in each formulation which can be set to zero. It 
is convenient to distinguish the coordinates constituting the vectors 
from the actual nodal degrees of freedom. The nodal degrees of freedom are 


denoted by q herein. The column vector containing all these freedoms, {q}, 

U 

may be directly identified with the vector {z} in equations 2 and 4. 


Symmetric Formulation 

For the symmetric analysis the airframe motions considered are pitch bend- 
ing of the fuselage, bending and' torsion of the wing and horizpntal tail, and 
fore and aft bending of the vertical tail. The aircraft was divided into 10 sub- 
structures having a total of l44 degrees of freedom. : The substructures employed 
and the corresponding degrees of freedom are identified in figure 19 and Table 7, 
respectively. 

Constraint equations .- Connection of the various substructures is achieved 
through the use of equations of constraint which mathematically enforce 
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deflection compatibility at the. interfaces of the substructures. The deflection 
vectors of adjacent substructures at their point of attachment must be equal 
when expressed in a common local coordinate system. The common systeias con- 
sxdered here are coincident vith the local coordinate systems of the contiguous 
substructures. The deflection vectors of adjacent substructures at their point 
of attachment are related by a diagonally partitioned rotational transformation 
matrix [rJ. Since small rotations can be treated as vectors the rotations at 
attachment points are related by the same matrix as the displacements. The two 
submatrices comprising [R] are thus identical and are denoted by [t]. For 
the general three-dimensional problem with six degrees of freedom associated 
with each station, [t] is of order 3 x 3 and [r] is of order 6 x 6. 

In the remainder of this section figure I 8 will be used in conjunction with 
auxiliary sketches, interspersed throughout the text , showing the substructure 
junctions and associated coordinate systems in order to aid in writing the 


equations of constraint. 



Sketch 1, Junction of fuselage beam #1 and fuselage beam #2. 
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The expression relating the equality of the deflections at the left end of 
beam H2 and the right end of beam #1 relative to coordinates local to beam ^2 
has the form 
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where the vectors {x}^ and {x}^ are identified with degrees of freedpm 
according to 
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and the coordinate rotation matrix [T^j 

r 

I cos 8^ 0 



sin 0^ 0 


is given by 


■sin 0^ 
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cos 0, 


(21) 
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Since there are no displacements or rotations out of the vertical plane of 
symmetry in a symmetric formulation y^, y^, a^, and are zero. 

The negative signs associated vith and in equations 20 have been 

introduced in order to have the usual definition of positive slope. This has 
heen done here merely for convenience. Several sign changes of this type will 
he Introduced during the course of this development for similar reasons. 
Expanding equation 19 using equations 20 and 21, the resultant constraint equa- 
tions at this junction are. given hy 

^2k '^11 ®1 - 9 ^ sin 0^ 

^12 ^11 ®1 ®1 ^ 22 ) 

“"^18 " "'^lO 

For brevity, the remaining constraint equations to be developed will not be 
written out piecemeal in the expanded form analogous to equation 22. A summary 
of all the constraint equations in expanded form will be given in an appendix. 



Sketch 2.- Junction of fuselage beam #2 and fuselage beam 
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At the junction of fuselage heams if 2 and §3 





Sketch 3*“ Geometric offset hetween fuselage and ving carzy-through. 


A rigid link is used to model the geometric offset between fuselage beam and 
the ving carzy-through beam-spring. The equations relating the deflections at 
the aft end of fuselage beam #2 to the deflections of the inboard end of the 
■wing carry-through beam-spring are then 
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Sketch 4.- Wing curry-through structure. 


Twisting of the wing carry-through structure was judged to be negligible in a 
syininetric mode of oscillation. Hence, the wing carry-through structure is 
taken to be rigid. This condition may be specified by the constraint 

'^50 - 1l)4 ° “ (=5! 

Extensional deformations of this member ar? also assumed to be negligible. 

Since = 0 (cf. eq^. 27) this implies the additional constraint 

= 0 (30) 
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fJketeh t).- Fuselage-wing jmiclipn. 


At the wing elastic azis/wing carry-through elastic axis junction 



where 



L|^ (Conversion axis) 


f wing/pylpn Junction. 
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The conversion axis is the centerline of the pylon conversion spindle which 
provides for tilting of the pylon. The geometric offset between the points of 
intersection of the wing elastic axis and the’ conversion axis with the wing tip 
rih (distance - .in the sketch'#6) is taken to “be rigid. The conversion axis, 
represented by the distance in sketch #6, .is also assumed tg be rigid. 

These rigidities are introduced by writing equations of constraint which relate 
coordinates to coordinates sequence of steps in establish- 

ing these equations .follow. 

An expression relating the coordinates at the last wing elastic axis 
station, intermediate (auxiliary) coordinates {x}^ is given by 



where {x}g and [T^^] are given by 





and 
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As pointed out earlier, the auxiliary coordinates {X} are not degrees of 
freedom but are introduced for convoni^jnee in arriving at the equations oi con 
straint. From sketch H6 above, the intermediate coordinate vectors {x}' 
and {X}^^ are seen to be related as 
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(37) 


If the conversion axis (length in the sketch //6) is taken to be flex- 

ible and treated as a beam-spring 

where {Q} is a column vector containing the degrees of freedom associated 
with the inboard end of the conversion axis beam-spring. The appropriate con 
straint equations would then have the matrix form 
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( 39 ) 


{X)g = [E,,] [L3] 


Similarly, at the outboard end of the, beam-spring representing the conversion 
axis 




(ho) 


vhere is a column vector containing the degrees oi’ fi-eedcm for the out- 
board end. If the conversion axis is taken as rigid {X}^^ and woiad 

not be identified vith degrees of freedom but would be auxiliary coordinates 
related according to 



Herein, the conversion axis is assi-uned to be rigid sc that equation 4l is 
applicable. 

•1 

% 

■ lit 

y 

Sketch 7 .- Junction of rotor shaft axis and e.niversioh axis. 



36 



The pylon structure, consisting of ths transmission/engine assemhly, is treated 
as rigid. The intermediate coordinate vectpr at the outboard end of 

the conversion axis is related to the coordinate vector {X}^^ situated on the 
proprotor shaft axis according to 
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The intermediate coordinate vector {x}^^ is related to the pylon center-of- 
gravity coordinate veqtor {X)^ by 
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•where 



On the basis of considerations in the preceding subsection and directly above 
{x}g and {X}^ are thus related accoij'ding to 


{X}g = [R^] [Lj] tEjl Wj 

The coordinate vectors -tx}g and (X}j_q are related as 



the beam-spring representing the mast, is given by 
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and {X}^ has already "been given in equjation i+5. 
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Sketch 8.- Geometry of pylon structure. 

The pylon structure, consisting of the transmission and an underslung engine, 
is supported by the conversion axis and the conversion actuator which attaches 
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to the forward end of the transmission case. This is assiuaed to be rigid. The 
flexibility of the mast, the portion of the drive shaft extending forward of 
the transmission case on which the proprotor is mounted, is represented as a 
beam-spring. However, the only elastic deformations of the mast which will be 
considered are vertical and lateral bending, motion in the other two directions 
being solely of the rigid body type. At the forward end of the mast: 
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To suppress the longitudinal and torsional deformations of the mast whlls allow- 
ing rigid-body motions in these directions write 


‘^101 - ^5 = ° 
^ 10 i + ■ '^8 ° 


( 50 ) 


These equations stipulate that the relative extensional and torsipnal deforma- 
tions between the ends of the mast are zero. 

The coordinates describing the proprotor disc, taken in the same 
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sense as { x ) 


Hence 




Assuming that the px'oprotor/huh comhination is rigidly fastened to the mast , the 
constraint equations are given by 



Sketch 9.- Attachment of vertical tail to fuselage "beam #3. 
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and 
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Since the geometric offset "between the fuselage and vertical tail elastic axes 
(distance L^) is taken to he rigid the intermediate coordinate vectors {x}^ 
and {x}'^ are related according to 
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and 



0 -cos9^ sin9^ 
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( 59 ) 


Assembling the series of transformation matrices indicated above the constraint 
equation is given hy 


{^>12 = [ R ^] { X }^3 


(60) 




At the jxinction of the vertical and horizontal tail elastic axes, 
{X}^' are related as 
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while and, {x}^^ are related according to 
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with {X}^ ^ and [T^] being given by 
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Egnations 6l and b4 ituply the niptrix relation 

{X}^ E [Rg] [ Ry ] {X}3_5 


(67) 
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This completes the derivation of the equations of constraint for the sym- 
metric formulation. For easy reference these equations are summarized in 
Table 8. Since there are k-9 constraint equations and l4i^ degrees of freedom, 
the matrix of coefficients of the constraint equations ([C] of eqs. 4 and 5) 
virill be of order 49 x 144. 

It should be noted that one of the constraint equations involving qj^^ is 
redundant. (Compare eqs. 8 and 10 in Table 8.) In practice, such redundancies 
cannot be avoided and inadvertently appear in the equations of constraint , 
resulting in equations which are not linearly independent. Since such redun- 
dancies are not usually identifiable by inspection, special consideration must 
■be accorded the constraint equations (ref. 22), The method advanced in refer- 
ence 19 obviates the need to treat such redundant equations of constraint in 
any special manner. Since this method is employed in the computational proced- 
ure implemented in the SUDM program which is used in the numerical example, the 
redundant equation will not (and need not) be deleted. 

pJ^o'Ee rties of beam components .- The mass and stiffness matrices 
for the beams em.ployed to represent the fuselage, wing, and empennage structures 
are established by substituting the data contained in Tables 1 to 3 into the 
appropriate matrix expressions developed earlier. The bending and torsion stiff- 
ness matrices follow from equations 6 and 11, respectively, while the lumped- 
mass inertia matrices, including wing and tail static unbalance, assume the form 
given in figure 9« Since the generation of these matrices using the data of 
Tables 1 to 3 xs stralghtfci-^-mra, the final numerical results are not shown. 

Structural proper cies of xigicl body components .- 

Wing c arry- through ; Although the elasticity of the wing carry-through 
structure is included, its inertial properties are treated as though it were 



rigid and combined with the inertia matrix of fuselage beam #2 by treating the 
center of gravity of the carry-through structure as rigidly attached to the last 
(right-hand) station of that beam. Assuming that the principal body axes at the 


center of gravity of the carry-through structure are parallel to axes x^, 

and the kinetic energy of the carry-through structure in terms of the 
3 

motion of its point of attachment is given by 


T 






( 68 ) 


where use has been made of the development shown in figure 10. Substituting 
equation 60 into Lagrange’s ecjuation and then making use of the appropriate data 
in Tables and 6 and equation 2U yields 
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as the matrix of additional terms to be added to the diagonal inertia matrix of 
fuselage beam ff2. 

Pylon; As indicated earlier, the rigid-body inertial properties of the 
pylon are defined relative to a coordinate axis system at the center-of-gravity. 
Hence, the pylon Inertia matrix follows directly from Table 4 and has the form 





where the minus sign associated with (cf, eq. 45) has been absorbed into 

the inertia matrix. The corresponding stiffness matrix is null. 

Proprotor. ihe propi’otor is treated as a rigid circular disc. Since the 
coordinate axes at the center of gravity of the proprotor disc are oriented 
such that they are principal axes the inertia matrix is diagonal and can be 
constructed directly from the data supplied in Table 4: 



(71) 


The matrix elements not shown are zero. The minus signs associated with IppQ 
*^111 51 ) *i 2 .VQ neon absorbed into the inertia matrix. The com- 

panion stiffness matrix is null since the proprotor hub is taken to be rigidly 
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attached to the mast. This matrix would not be null if the hub were spring- 
connected to the mast j as will be demonstrated below. 

Structural properties of spring component s .- 

Wing carry-through: Since x^ has been set to zero (cf. eg. 27) the beam- 

spring stiffness matrix is of order 11 x n and is given by that in figure 5 
with the first row and column deleted. This matrix is shown in figure 20 along 
with the ordering of the degrees of freedom. Tables 5 and 6 contain the neces- 
sary data to evaluate the terms of the matrix. Constraint equations 29 and 3G 


in conjunction with the 11 x ii stiffness matrix given in figure 20 imply that 
twisting and both the elastic and rigid-body axial motions of the beam-spring 
representing the wing carry-through structure are removed as degrees of freedom, 
in accordance with earlier discussions. 

Pylon mast: The stiffness matrix of the pylon mast beam-spring is given 

in figure 21 along with the ordering of the degrees of freedom. Again, Tables 5 
and 6 contain the data required to evaluate the individual terms of the matrix. 

Some comments are included here to indicate the manner of treating a 
proprotor/hub assembly which is connected to the mast by springs. For illustra- 
tive purposes assume that the proprotor/hub combination is allowed to flap 
longitudinally and laterally with respect to the mast, the flapping motion being 

restrained by rotational springs K ' and IC , respectively. Since longi- 

^1 ^1 

tudinal and lateral angular motions of the proprotor disc relative to the mast 
are permitted constraint equations 39 and 4o in Table 8 must be deleted. Recall- 
ing the treatment of springs in the spring-mass analogy employed earlier to 
illustrate the manner of including fuel slosh in a launch vehicle vibration 


analysis, an expression for the strain energy stored in the springs is required. 
The appropriate expression here is given by 

50 




(72) 



as the matrix of spring terms to he added to t^ie partitioned system stiffness 
matrix, [K], For the coordinate pumhering shown in fi^iiire 19, it is to he 
noted that the off-diagonal terms in equation 7^ will couple the stiffness 
matrices for the pylon mast and proprotor suhstruqtures, 

Anti-Symmetric Formulation 

The considerations related to the symmetric formulation have illustrated 
the manner of establishing substructure mass and stiffness matrices and the 
mechanics of setting down equations of constraint, A corresponding derivation 
for the anti-symmetric case would, for the mpst part, he repetitious, For this 
reason a summary-type treatment , listing only final results which are different 
from the symmetric case, is giyen here. 
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For the anti-symmetric analysis the airframe motions considered are side 
bending and torsion of the fuselage and vertical tall and bending and torsion 
of the wing and horizontal tail. The same substructures and stations used in 
the symmetric analysis are employed here, leading to a problem having I 65 degrees 
of freedom. The substruct\ares are ‘identified in figure 22 and the pertinent 
degrees of freedom in Table 9. 

Constraint equations .- For anti-symmetric motions of the airframe the twist 
of the wing carry-through structure was judged to be important and the beam- 
spring representing it now also being permitted freedom in twist. Although 
ertensional deformations of the wing carry-through are still negligible the 
3-bility to translate axially as rigid body must be provided for. The final 
constraint equations, 48 in number, are summarized in Table 10. The matrix of 
coefficients of the constraint equations ([C] of eqs. 4 and 5) is thus of 
order 48 x 165. 

Structural properties of beam components .- As in the symmetric case the 
mass and stiffness matrices for the beams employed to represent the fuselage, 
wing, and empennage structures are established by substituting the data in 
Tables 1 to 3 into the appropriate' matrix expressions developed earlier. Again, 
because the generation of these matrices is straightforward the final numerical 
results are not shown. 

Struct\nral properties of rigid body component s .- 

Wing carry-through: Here the vector is given by 
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Pylon: The pylon inertia matrix is identical to that given by equation 70 

for the symmetric case, with the degrees of freedom shown there, thro^gh 

^ 94 » replaced by through ^espeptiyejLy. 

Proprotor: The proprotor inertia matrix is identical to that giveq by 

equation 71 for the symmetric case, with the degrees of freedom shown there, 
qj _07 through ^2.12’ ^^Pl^’Ced by q^^j^ through respectively. 



Wing carry-through: Since axial rigid-body motion of the parry-through is 
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included in the anti-symmetric formiaation the heam-spring representing it has 
a stiffness matrix of order 12 x 12. Its form is the same as that given in 
figure 21. However, the degrees of freedom shown there, through 1^06 ’ 

are replaced hy through respectively. Tables 5 s-nd 6 contain the 

data needed to evaluate the individual terms of the matrix. 

Pylon mast: The pylon mast beam-spring stiffness matrix is also given by 

figure 21 if through q^^g are replacecj by q ^^0 through qQ_ 23 » respec- 

tively. Again Tables 5 and 6 provide the daba required to evaluate the indi- 
vidual terms of the matrix. 

Niimerical Resu;^.ts 

The direct method of analysis p-s embodied in the SUDAN program was used to 
calculate the symmetric and anti—symmetric free— free modes and frequencies for 
several pylon tilt angles. A STommary of these frequencies for the first five 
elastic modes is given in Table 11. The first five symmetric modes for the case 
in which the pylons are tilted fully forward ~ 90^) are sketched in fig- 
ure 23. For reasons of both pictorial clarity ai>d convenience, the changes in 
slope of the fnselage elastic axis and the offsets of the wfng and vertical tall 
elastic axes from the fuselage elastic axis are not shown in the sketches, Ifc 
should also be noted that while all degrees of freedom contribute to a coupled 
mode shape the shape is oftentimes dictated by the motion in a relatively few 
of the degrees of freedom. Such is the case here and only the predominant 
motion is indicated in the mode shape sketches given in figure 23. 
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CONCLUSIONS 

The manner of representing a flight vehicle structure as an assembly of 
beam, spring, and rigid-body components for vibration analysis has been des- 
cribed, The development was couched in terms pf a substructures methodology 
which is based on the finite-element stiffness method. Basic concepts were 
first introduced through several specific qualitative examples which ipcluded 
the modeling of wing structures, external stores, pylons supporting engines or 
external stores, and sprung masses associated witji launch vehicle fuel slosh. 

A detailed numerical example consisting of a tilt-rotor VTOL aircraft was also 
given to provide an illustration of the procedures for modeling a structure, the 
manner of forming the substructure mass and stiffness matrices, and the mechanics 
of writing the equations of constraint which enforce deflection compatibility at 
the Junctions of the substructures* Since many structures, or selected compo- 
nents, can be represented in this manner for vibration analysis, the modeling 
concepts described and their applications in the numerical example shown jshould 
prove generally useful to the dynamicist. 
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TABLE 1.- FUSELAGE DISCRETIZATION 


Station 

Local 

Coordinate 

Position 

Mass 

Torsional 

Inertia 

El 

Vertical 

.X 

Side 

GJ 

(cm) 

(kg) 

(kg-cm^) 

(kB-cm^) 

(kN-cm^) 

(kN-cm^) 

Beam #1 

, 



. 



1 

0.0 

18.2 

5.65 X 10^ 

9.75 X lo"^ 

8.61 X 

10.33 X 10*^ 

2 

88.9 

21.0 

11.30 

22.96 

21.52 

24.11 

3 

152.4 

52, 5 

45.19 

45.92 

43.05 

57.40 

k 

215.9 

252.2 

282.46 

91.83 

80.35 

114.79 

5 

304.8 

378.3 

175.12 




Beam #2 







1 

0.0 

378.3 

175.12 X 10^ 

229.58 X 

129.14 X 

200.89 X 10 

2 

101.6 

924.7 

420.30 

390.29 

218.10 

246.80 

3 

190,5 

357.2 

271.16 

459.17 

269.76 

252.54 

k 

254.0 

42.0 

124.28 

476.39 

286.98 

252.54 

5 

304,-8 

168.1 

259.86 

482.13 

298.46 

254.54 

6 

355.6 

147.1 : 

175.12 




Beam #3 




1 , 

0,0 

147.1 

175.12 X 10^ 

482.13 -x lo"^ 

299.89 X lo”^ 

252.54 X lo"^ - 

2 

86.9 

328.9 

412.39 

476.39 

298.46 

252. 54 

3 

l60,0 

507.9 

429^34 

418.99 

252.54 

252.54 

4 

266.7 

1849.3 

463.24 

258.28 

172.19 

232. 45 

5 

381,0 

447.6 

463. 2l 

129.14 

114.79 

172.19 

6 

477.5 

126.1 

101.68 

86.09 

80.35 

114.79 

7 

596.9 

61.3 

67.79 

51.66 

51.08 

74.61 

e . 

729.0 

105.1 

45.19 

— 




I Ictal Mass: Beam #1 = 722.2 kg, Beam #2 = 2017. U kg. Beam #3 = 3573.3 kg 




TABLE 2.- WING DISCRETIZATION 


Station 

Local 

Coordinate 

Position 

Mass 

Static 

Unbalance* 

Torsional 
Inertia 
About e.a. 

El 

Vertical 

(klJ-cm^) 

El 

Chord 

GJ 

(cm) 

(kg) 

(kg-cm) 

(kg-cm2) 

(kN-cm^) 

(kN-cm^) 

1 

0.0 

154.3 

-2304. 

4.4o X 10 ^ 

4.30 X 10 ® 

10.04 X 10 ^ 

3.53 X 10® 

2 

86.4 

123.5 

-2882. 

4.68 

3.56 , 

8.60 

2.98 

3 

190.5 

90.7 

- 2188 . 

2.78 

2.64 

6.89 

2.35 

h 

317.5 

113.5 

- 2478 . 

2.61 

1.72 

5.16 

1.72 

5 

439.4 

62.5 

- 2882 . 

2.84 

1 . 32' 

4.30 

1.32 

6 

530.9 

259.2 

-4i46. 

5.94 

1.19 

4.08 

1.20 

7 

563.9 

0.0 

0 . 

0.00 

— 

— 

— 


Total Mass = 803.6 kg 


^Negative static unbalance indicates that section center of gravity is aft of section 
elastic axis. 


TABLE 3 .- EMPEMAGE DISCRETIZATION 



LCtal Mass = 


C.O 
92 . ^ 
181.6 

274.3 

363.3 


94.9 kg 


Mass 

Static 

Unbalance* 

Torsional 
Inertia 
About e . a . 

El 

Side 

El 

Chord 

(kg) 

I 

(kg-cm) 

( kg-ciri^ } 

(ki^-CE‘^) 

(klJ-cm^) 

15.1 

22.9 

21.8 

20.4 

i4.6 

751.7 
! 622.3 

564.9 
467.1 
353.6 

1.58 X 10 " 
1.23 
.98 
.72 
.49 

5.02 X lo'^ 
1.86 
.57 
.29 

4.36 10 ® 

2.30 

.98 

.52 



3. 73 X IC' 

2.27 

1.15 

.52 


Horizontal 

Tail 

1 

0.0 

0.0 

0.0 

0.0 

1.29 X 10^ 

1 

.86 X 10 ® 

1.15 X 10^ 

2 

88.9 

0.0 

0.0 

0.0 ^ 

.57 

.43 

.79 

3 

k 

165.1 

63.0 

0.0 

.25 X 10 '' 

.34 

.21 

.29 

254.0 

0,0 

0.0 

0.0 

— 

— 

— 

Total Mas: 

3 = 63.0 kg 






1 


c 

>T! fe 

o t 
o 

gd W 
s 


^Positive static unbalance indicates that section canter of gravity is forward of sectioi 
elastic axis. 




























TABLE U.- INERTIAL PROPERTIES OF 

COMPONENTS TREATED AS RIGID 

PROPROTOR 


Mass 

637.4 Fg 

Flapping Inertia 


Longitudinal 

1.6087 10^^ kg-cm^ 

Lat eral 

1.6087 lo"^ kg-cm^ 

Polar Inertia 

3.2174 X io'<' kg-cm^ 

PYLON 


Mass 

1269.7 kg 

C. G. Inertia 


Pitch 

7.8524 X 10^ kg-cm^ 

Roll 

4.2934 X 10^ kg^cm^ 

Yaw 

6.9146 X 10^ kg-cm^ 

WING CARRY-THROUGH * 


Mass 

154.1 kg 

C. G. Inertia 


Pitch 

8.1936 X 10^ kg-cm^ 

Roll 

Not avail a.ble 

Yaw 

Not available 

* 

Treated as rigid inertially but not elastically 

6l 




'TABLE 5.- STIFFNESS PEOPERTIES OP COMPONENTS TEEATED AS SPEINGS 
Wing Carry-Through* 


El 


vertical 


El. 


lateral 


GJ 


= 4.5056 X 10 ^ NN-cm^ 
= 10.5898 X 10® NN-cm^ 
= 3.6733 X 10 ® kN-cm^ 


Pylon Mast 


El 


vertical 


El 


lateral 


7.1745 X 10’^ kN-cm^ 
7.1745 X lo"^ kN-cm^ 


*AE is taken to be infinite 

*(• 

Both AE and GJ are taken to be infinite 


62 



TABLE 6.- GEOMETRIC QUANTITIES EMPLOYED IN VIBRATION ANALYSIS 


ANGLES 


0 


1 



0 


4 


0 

0 


5 

6 


LENGTHS 


■^1 



9 » 5 ° (.1658 Radians) 
5.6° (. 0977 ^^ Radians) 
4.5° (.07854 Radians) 
Varied 

25° (.4363 Radians) 
15° (.2618 Radians) 


117.1 cm 

106.7 cm 

43.2 cm 

43.2 cm 

79.2 cm 

48.3 cm 

60.4 cm 

71.1 cm 

38.1 cm 
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TABLE 7.- IDENTIFICATION OF DEGREES OP FREEDOM POE 


SYMMETRIC VIBRATION ANALYSIS 


FUSELAGE 


Beam #1 

- S 

Displs. \ 

% ~ ^10 

f 

Slopes ) 

‘^11 

Axial rigid body 

Beam §2 

^12 “ '^l? 

Displs. ) 

^18 ” ^3 

Slopes / 

^2h 

Axial rigid body 

Beam ^3 

^23 ~ ^32 

Displs. 1 

*^33 " ‘^UO 

Slopes j 

'I41 

Axial rigid body 

Wing Carry-Through Structure 

142 “ ^^2 

WING 

S 3 S 9 

Displs. 1 

So S 6 

Slopes ) 

St " Ss 

Displs. 1 

S 4 “ Sc 

Slopes ) 

Si " St 

Torsion 

S 8 

Axial rigid body 


Vertical Bending 


Vertical Bending 


Vertical Bending 


Vertical Bending 


Fore and Aft Bending 


6)4 



TABLE T*“ Concluded. 


Assembly 

Rigid Body Translations and Rotations 


Rigid Body Translations and Rotations 


PYLON 

Transmission/Engine 
^89 ” 

Mast 

^5 “ ^^106 

PROPROTOR 

*^107 “ ^112 

EMPENNAGE 
Vertical Tail 

“^iis '^117 
'^118 " '^122 
%23 

Horizontal Tail 

%24 " ^127 
^128 " '^ISl 
^^132 " '^135 
%36 “ ^^139 
‘^li^O “ ‘^lit3 


Bispls I 
Slopes * 

Axial rigid body 

Displs . 

Slopes 
Displs . 

Slopes 
Torsion 
Axial rigid body 


Fore and Aft Bending 


Vertical Bending 


Fore and Aft Bending 
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TABLE 8.- CONSTRAINT EQUATIONS FOR SYM4ETRIC VIBRATION ANALYSIS 


1. 

hi 

cos 9^ 

- I5 

sin 

01 - 

hk ~ 

0 

2. 

“^ll 

sin 9^ 


cos 

01- 

' h2 ~ 

0 

3. 

"‘‘^10 ^ '^iS 

= 0 





L. 

%1 

cos 0p 

" ^^25 

sin 


'■^24 

^ 0 

5. 

%1 

sin 

Cl 


cos 

e. 

” -17 "" 

■ 0 

6» 

“^33 ^^23 

= 0 





T. 

942 

+ Li q^. 

1, + <hk ■' 

0 



8. 

-^1 

Qi ,_ = 0 

^45 





9. 

%3 

" 9it - 

' 0 





10. 

“94 c 

. = 0 
) 






11. 

944 


^ 0 





12. 

%6 

= 0 






13. 

^^50 

- 944 = 

' .0 





14. 

947 

= 0 






15. 

‘^88 

cos 9^ 

“ *^67 

sin 

S 

- 947 = 

0 

16. 

^^88 

sin 9 2 

hi 

cos 


%8 “ 

0 

IT. 

Ss 

- 944 = 

^ 0 





18. 

% 

cos 9^ 
3 

^ *^60 

sin 

®3 

wj rj xzi 

HiCJ p. 

0 

19. 

'ISl 

sin 6^ 

ho 

cos 

®3 

~ ^51 == 

0 

20. 

9j4 

^52 

0 





21. 

'^89 

cos 9^ 
3 

Sto 

sin 

9., 

CQ3 6^ • 

” ^^91 


- 

Lg sin 

6^ CCS 9i 

»v *" 

]- 

9o -> 

9:i 1 p 

cos 


+ 


cos 6,_, - 

cos 

^4 (H ■ 

sin 9 


3^ sin 0 ^ + sin sin 

.3 - sin 0 ^(l 


e 


'4 


sin 9^ + cos 6 ^ 


3 . 

cos - qgg = 0 
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TABLE 8.- Continued 


22. 


23. 

24, 


sin + q^Q cos 6^ cos 6j^ - cos 0^ 3in 9^ + ^^92 [^5 ®4 


Lg cos 6^ cos 


*^3 


•L/- sin 6^ + sin 6i. (L,, sin 0^ - Li. cos 


41^3 


q^]^ sin 9^ + cos ®3 “ ^4 ®3 


$)_ “ ^^73 ° 


%2 ®3 “ '^c.3 ®3 ■ ^91^ ®3 ®4 " ‘ 3-87 ° 

.] 


<150 ®i, * V ®i, + 152 



— 

^93 ^4 

cos 0j^ - 

I9I, H 01, - 

25. 

I92 ®3 * 

q^3 cos 

®3 cos 

26. 

"*^93 

sin 0j^ 

+ q^;^ cos 

0 

11 

0 

CO 

r 

CD 

CM 

"'^96 

^'7 '^lo'o ^89 

+ Lg q^g = 0 

28. 

^^95 “ 

■ <190 ^ 

^6 *^99 

0 

29. 

I97 - 

■ ^ S9 

" '^91 " 

0 

30. 

^99 - 

^92 

0 


31. 

“'^98 

+ ^93 = 

0 


32. 

^100 

- S4 = 

0 


33. 

‘^101 

- I95 = 

0 


34. 

^^104 

■ Sq " 

0 


35. 

^^107 

" '^101 

= 0 


36. 

‘^108 

^102 

= 0 


37. 

'^109 

" %03 

= 0 


38. 

“'^lic 

1 ^io4 

= 0 


39. 

“'^111 

. "" ^105 

= 0 


4o. 

^112 

■ '^106 

= 0 


4l. 

^113! 

-cos 0g 

cos 0f- + 

5 

■ sin Sg sin 0^ 



sin 0g 

cos 0^J + 

•Ilia S ®a 


L^ - L^ cos 9)^ 


J-59 


Lg sin 61,1 


^123 


cos 9g sin 9^ 
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TABLE 8.“ Concluded 


k2. [sin cos 0^ + sin 0^ cos OgJ + -sin 0^ sin 0^ 

+ co= 02 =os 0jj - Lg sin 0^ - qgg = 0 

its. + q^g = 0 

ii’4. cos 0g + sin Gg = 0 

^5- -q.i4i^ cos 0^ sin 0g + q^22 ^5 ^6 ®5 " I115 = ° 

U6. sin 0^ sin Og - sin 0^ cos 0g + cos 0^ - = 0 

kl. cos 0g - q^oRsin 0g + --- 0 

^8* -^il^o 8^ sin 9g - q^^g cos 0^ cos 9g -!■ q^^g sin 0^ = 0 
^l4o win Gg 4" qj^.-)Q ^in 0^ ccs 6g + 6^ — 0 
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TABLE 9.- IDENTIFICATION OF DEGREES OF FREEDOM FOR 
ANTI-SYMMETRIC VIBRATION ANALYSIS 


FUSELAGE 

Beam 

\ - S 

% " 

“ ^15 

Beam 

'^16 “ ^21 

922 “ ^27 

q-28 - 933 
Beam ^3 

^34 " ^4l 

%2 ~ ‘^49 

'I 50 ■ St 

Wing Carry- Through 
SS “ %9 

WING 

50 “ S6 
St " S 3 
S4 “ So 

51 “ St 
Ss “ So4 

S 05 


Displa. 

Slopes 

Torsion 

Displs. 

Slopes 

Torsion 


Displs. 

Slopes 
Displs . 

Slopes 
Torsion 
Axial rigid body 


Displs 

Slopes 

Torsion 


1 


Structure 




Side Bending 


Side Bending 


Side Bending 


Vertical Bendipg 


Fore and Aft Bending 
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TABLE 9*“ Concluded 


PYLON 

TransmisEion/Engiiie Assembly 

" '^121 Rigid Body Translations and Rotations 

Mast 

^112 " %23 
PRQPROTOR 

^12h ~ ^129 Rigid Body Translations and Rotations 

EMPEMAGE 
Vertical TaiX 


'^130 ~ Displs. 

^^135 ~ '^139 Blopes 

%40 " \kh • Torsion 

Horizontal Tail 

^1^5 “ %U8 Displs. 

'll 1+9 ~ %52 Slopes 

'^153 “ ^156 Displs. 

'^157 " ^l6o 



^l6l “ '^161^ 


Torsion 


Axial rigid body 


Side Bending 


Vertical Bending 


Fore and Aft Bending 




TABLE 10.- CONSTRAINT EQUATIONS FOR ANTI-SYMMETRIC VIBRATION ANALYSIS 


2 . cos 0 j_ - q^Q sin \ ^ 0 

3 . q^^ sin 0 ^ + q^^ cos 0 ^ - q^^ =? 0 

131 . - 121 “ ° 

5. qjQ cos Sg - sin 6g - ^33 = 0 

6 . q^(j sin Sg * q^g cos Sg - qg^ = 0 

T. qj5 + qg^ = 0 

8. qjg - 1-3 qgg - qgi = 0 

. s- IfiO = 0 

10 . qgg + q^^ = 0 

11 . qg^ =0 

12. qg 3 - qg.^ = 0 

ifiit - I58 “ ° 

111. q^Qj cos 63 - qg^ sin 83 - qg4 = 0 

15. qj^Qj sin 63 + qg^ cos 83 - qgj = 0 

’to - 166 = ° 

IT. q^Q cos ijrj sin 0^ - qg^ = 0 

18 . q^g sin 0 ^ - cos 83 “ agg = 0 

19 . q^^ - qg^ =0 

20 . q^Qg cos 0 g + q^Q^ sin 0 g cos 0 ^^ - q^^g sin 0 g sin 9 |^ 

■" ‘‘logo's ®3 \ - *-6 ®3 \1 ■ »llo[’^6 ®3 

" ®3 * ^3 ^ 3 )^ ■*" ®3 " 83 

- Lg cos 83 )] - q^05 = “ 
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TABLE 10.- Continued 


21. 

-‘I106 S ^ *^107 ®3 % 

- qio8 83 sin 0J^ 


* ^109 [^5 ®4 ®3 ” ^6 ®3 ®4j ^ 9-110 "^6 ®3 


+ sin 011^^3 sin 6^ - cos 

J '*’ *^111 L 5 ®3 ^4(^3 ®in ®3 


- cos 03)j - q^o = 0 


22. 

qi09 cos 03 - q^^Q sin 83 cos 0j^ - 

=^'' ®3 ®4 - %oit “ ° 

23. 

I10Y sin ej, -^ q^gg cos -. q^Qjji 

13 - Lg 00s 6,^ - Lg sin ej 


- 9iio Lj, cos - qj^^^ Lj, sin 

1 ^ 6 ‘° 

2 k. 

-q^09 sin 03 " ^hq cos 83 cos 0j^ 

- qiii cos 83 sin 6 ^ + 1Q3 = 0 

25. 

-9iio sin 0j^ + q^^^ cos 0^^ - = 

0 

26. 

"'^113 ^7 ^^117 “ '^106 '^ H ^^115 

0 

27. 

^112 ‘ ^107 ^6 ^116 = ° 


28. 

%l4 " ^7 ‘^116 “ ^108 ° 


29. 

‘^116 “ ‘^109 “ ° 


30. 

-‘^115 ^110 = ° 


31. 

^^117 ~ ^111 ° 


32. 

%18 " 9-112 == ^ 


33. 

^121 *■ '^115 ^ ° 


34. 

%24 " ^118 "" ° 


35. 

%25 " ^^114 ~ ° 


36. 

^^126 " ‘^120 ° 



37 . - qj _27 ■^ '^121 ^ ° 

"'^128 ‘^122 ^ ° 

39 . q^29 “ 9-123 “ ° 

^0- 9.13Q - ^ 13 ^ cos 6^ + sin 0^ - qj^^ = 0 


72 



TABLE 10,- Concluded 

^1* ^IS^Ijcos 02 cos 0^ + sin 0g sin oj + jcos 02 sin 0^ 

+ sin 02 cos 0^j - = 0 

^ 2 . ^2 ®5 “*■ ®2 ^ 5 ] '*’ ‘^ 14 q ( t ®^’‘^ ®2 ®5 

+ cos 02 cos 0^j - = 0 

43. q^g5 cos 0g + q^53 sin 0^ - q^^32 = 0 

44. -q^g^ cos 0^ sin 0^ + q^^3 cos 65 cos 0g + q^j^^ Sin 6^ * 0 

45. q^g^ sin 0^ sin 0^’ - q^^3 sin 0^ pos 0g + cos 6^ = 0 

46. q^g^ cos 0g - q^i^^ sin 0^ = 0 

"^161 ®5 ®6 ” '^149 ®5 ®6 ®5 " ‘^IST 

48. q^g^ sin 0^ sin 0^ + q^j^^ sin 0^ cos 0g + q^^^ cos 0^ - q^^g 
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TABLE 11.- FREE-FREE ELASTIC MODE FREQUENCIES OF MODEL 266 TILT-ROTOR (Hz) 


Symmetric Mode Frequencies 


\ 0 
ModV 

0 ° 

15° 

30 ° 

45° 

0 

0 

■ 75 ° 

90° 

1 

2.114 

2.138 

2.168 

2.199 

2.224 

2.235 

2.229 

2 

3.473 

3.488 

3.487 

3.470 

3.304 

3.399 

3.360 

3 

5.727 

5.698 

5.661 

5.628 

5.610 

5.617 

5.644 

4 

7.235 

7.378 

7.646 

8.021 

8.694 

8.895 

9.016 

5 

11.490 

11.509 

11.522 

11.523 

10.802 

11 . 046 

10.661 


Anti-symmetric Mode Frequencies 


\ 9 

0 ° 

15 ° 

— 

30 ° 

45° 

60° 

75 ° 

90° 

1 

3.882 

3.882 

3.884 

3.883 

3.875 

3.856 

3.825 

2 

4.602 

4.699 

4.801 

4.891 

4.936 

4.906 

4.818 

3 

5.702 

5.645 

5.569 

5.492 

5.446 

5.464 

5.528 

4 

7.256 

7.347 

7.373 

7.339 

7.294 

7.251 

7.214 

5 

7.633 

7.624 

7.775 

8. 073 

8,456 

8.860 

9.136 





















Beam element 1 
EI^,GJ^.AEj, ,L^ 


Beam element 2 
Elg.GJg.AEg.Lg 


Figure 3.- A two-element Beam substructure. 



Figure k.- Sign convention for end deflections of massless uniform 

beam segment (beam-spring). 
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Figure 6. 


opi xiig 

Rigid Body 


- "stick" model of transport aircraft for symmetric 
vibration analysis. 


BSJRODUCBIUH®®' 
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I Substructure #1 - forward fuselage beam 


' Substructure §2 - aft fuselage beam 


Axial rigid body 



/ Substructure #4 - wing beam 


Substructure |5 - vertical tail 


Substructure #6 -fiorirontal tail 


Axial rigid bcdy- 


Substructure *8 - inboard nacelle ^ 


Substructure #7 - 
Inboard pylon 


Substructure *9 - oulboanJ pylon 


Substructure #10 - outboard nacelle 




Block diagonal conpc sition of matrices [M] and [k] for the -beam model of Figure 6. 



Figure 8.- Identification of mass coupling terms arising from wing static tmtalance. 



Figure 9.~ Coupled bending-torsion mass matrix corresponding to the wing of figure 8 






Figure 10.- Manner of accommodating the rigid-body inertial properties 
of a center-line fuel '.anlc into the inertial matrix of the fuselage 
beam. 
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CX3 

-c- 


liL, 




Figure 11.- Inertia matrix of the fuselage benm/fuel 









(b) Springs having both ends free to move. 


Figure 12.- Some posible beam/spring arrangements which can arise from 

a structural idealization. 
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(a) Springs ’ having one end tied to ground. 


Figure 13.- Form of stiffness matrices corresponding to the spring components 
in several of the beam/ spring arrangements of figure 12. 
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Figure l4. Use of tlie spring-mass analogy to simulate the dynamics of sloshing 

propellant in a natural “mode anaJ.ysis. 








(a) Side view 


Figure 18. Coordinate systems used to write the equations of deflection compatibility 

between adjacent substructures. 




Figure 18 .- Concluded 
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Figiire 19 •- Block, composition of matrices [mJ and [K] for symmetric analysis of 'beam 

model shown in figure IT. 
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Figure 22.- Block diagonal composition of [M] and Tk] for anti-symmetric analysis of 

*KcicjTn QinnVn *1 Tl IT 





